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We analyze the properties of four-quark condensates and scalar susceptibilities in the meson gas, 
within finite temperature Chiral Perturbation Theory (ChPT). The breaking of the factorization 
hypothesis does not allow for a finite four-quark condensate and its use as an order parameter, except 
in the chiral limit. This is rigorously obtained within ChPT and is therefore a model-independent 
result. Factorization only holds formally in the large A^c limit and breaks up at finite temperature 
even in the chiral limit. Nevertheless, the factorization breaking terms are precisely those needed 
to yield a finite scalar susceptibility, deeply connected to chiral symmetry restoration. Actually, we 
provide the full result for the S'(7(3) quark condensate to NNLO in ChPT, thus extending previous 
results to include kaon and eta interactions. This allows to check the effect of those corrections 
compared to previous approaches and the uncertainties due to low-energy constants. We provide a 
detailed analysis of scalar susceptibilities in the SU (3) meson gas, including a comparison between 
the pure ChPT approach and the virial expansion, where the unitarization of pion scattering is 
crucial to achieve a more reliable prediction. Through the analysis of the interactions within this 
approach, we have found that the role of the a resonance is largely canceled with the scalar isospin 
two channel interaction, leaving the piJlQ) as the main contribution. Special attention is paid to 
the evolution towards chiral restoration, as well as to the comparison with recent lattice analysis. 
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I. INTRODUCTION 



Chiral symmetry restoration T] is a very relevant ingredient in our present understanding of hadronic physics under 
extreme conditions of temperature and density and has been one of the main motivations for the development of the 
heavy-ion and nuclear matter experimental programs, which are still producing new results in facilities such as RHIC, 
CERN (ALICE) and FAIR. In parallel, the improvement of lattice data at finite temperature, performed by different 
groups [SHI], has contributed considerably to clarify the main properties of the chiral symmetry transition, which is 
believed to take place in the same range as the deconfinement one. Nowadays, there is a fair consistency between 
lattice simulations performed with different methods, pointing towards a crossover-like transition for Nf = 3 (2-1-1 
flavors in the physical case), becoming of second order for Nf = 2 (in the 0(4) universality class) and first order in 
the degenerate case of three equal flavors. This behavior corresponds to vanishing baryon chemical potential and the 
transition temperature lies within the range ~ 150-175 MeV. It is important to remark that in the physical 2-1-1 
case analyzed on the lattice, the transition being a crossover means that one should really talk about a transition range 
rather than a critical temperature, and that range can be established by looking at different order-like parameters, 
which can give different values for T^. The chief parameters used in lattice analysis are the quark condensate and 
susceptibilities, deflned as first derivatives of the quark condensates with respect to quark masses. Susceptibilities 
measure fluctuations of the associated order parameter and can be expressed in terms of current correlators. Thus, 
the scalar susceptibility, related to the quark condensate, is expected to grow faster just below the transition. 

It is important to provide an accurate analytical description of the physics below the chiral transition, to com- 
pare with experimental data and to confront the lattice results in the continuum. On the one hand, a particularly 
useful approach has been the Hadron Resonance Gas (HRG), which has proven to be quite successful to describe 
thermodynamic quantities, i.e., those that can be derived directly from the free energy density, as compared to lattice 
data [8Hl0i|. Within the usual HRG approach, the free contribution of all known physical hadron states to the parti- 
tion function is considered. Models including resonance widths and hadron interactions improve the HRG approach 
description of hadron production experimental data [TTJ [H] and lattice results TS] . 

On the other hand, Chiral Perturbation Theory (ChPT) [I4HI6] allows to describe the low-temperature meson gas 
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[T7] in a model-independent and systematic way for SU{Nf)L x SU{Nf)ii — SU{Nf)v symmetry breaking with 
Nf = 2,3 light flavors. In particular, interactions can be included in the most general way compatible both with the 
underlying symmetries and with meson-meson scattering data. Although ChPT includes only the lightest degrees 
of freedom {tTjK,!]) and hence is limited to low and moderate temperatures, it provides model-independent results. 
In fact, in early studies of the partition function ,17J , extrapolations of the condensate and other quantities to the 
transition region give a qualitatively reasonable description of the relevant physics. In addition, there have been 
some relevant developments on the use of ChPT results in the description of the meson gas properties. For instance, 
it is known that the Inverse Amplitude Method P^l - BU] . which is a dispersive method to unitarize ChPT without 
introducing spurious parameters, allows the generation of light meson resonances ( the a or /o(500), the p{770), ...). 
Within that unitarized ChPT scheme, useful results have been developed for the temperature and density dependence 
of the lightest resonances and their connection to chiral restoration and chiral partner degeneration as well 

as a phenomenologically successful description of transport coefficients 24J. It is particularly relevant to recall the 
virial expansion [fTl l25H30j which allows to parametrize efficiently the effect of meson interactions in the partition 
function within a dilute gas description, valid below the transition. States of more energy are weighted by Boltzmann 
factors and then become more relevant as the system approaches the transition. Thus, in this approach, it is more 
important to include accurately the interaction of the lightest mesons, e.g. via unitarization, while the heavier ones 
can be added as free states. Actually, the HRG approach with just free states is nothing but the leading order in the 
virial expansion. An alternative approach, also within unitarized ChPT, is to take into account the temperature or 
density dependence of the phase shifts, which would make the a resonance pole move towards the real axis precisely 
as a signature of chiral restoration [531 123] • 

In this work we will explore some additional properties of the meson gas within the ChPT framework, regarding in 
particular quark condensates and susceptibilities. Our study will be deeply connected to the analysis of four-quark 
correlators of the type {T{qq){x){qq){0)) and the factorization hypothesis, thus extending our previous work at zero 
temperature [3T]. This hypothesis states that four-quark condensates factorize into two-quark condensates squared 
{{qq)^) ^ (99)^ with the same quantum numbers. We have shown in [31] that this hypothesis fails to next to next to 
leading order (NNLO) in ChPT. Here, we will show that the same holds at finite temperature, preventing the use of 
four-quark condensates as order parameters since the factorization breaking terms diverge. In the derivation we will in 
turn obtain the SU{3) thermal quark condensate in ChPT including all the relevant meson interactions, which extends 
previous calculations of the condensate to this order which considered free kaons and etas |17| . The non- factorizing 
scalar four-quark correlator gives rise to the scalar susceptibility, allowing for a direct check of the calculation. Once 
the connection with factorization is established, we will perform a detailed analysis of the scalar susceptibilities in 
the virial approach, with and without including unitarized amplitudes, extending previous works in the literature and 
serving as a test of the robustness of the ChPT results, which together with their model independence makes them a 
useful prediction for low and moderate temperatures below the transition. In this respect, it is particularly relevant 
to note that unitarized ChPT is able to provide a relatively good description of the quark mass dependence of the 
resonant states obtained through unitarization [321 - I34] . particularly robust in the case of the lightest scalar. 



The paper is organized as follows. After fixing our notation and definitions (section II A), in section II B we present 



our calculation of the relevant four-quark correlators for two and three flavors at flnite temperature. The details are 



given for Nf = 3. The factorization hypothesis is then examined in section II C[ where we also comment on the 
large- iVc limit. In section IIP we establish the connection with the scalar susceptibility. Different effects in the ChPT 



condensates and susceptibilities are discussed in section |III[ while section |IV| is devoted to the analysis within the 
virial approach, to compare with previous analysis and to study the role of interactions. We pay special attention to 
the comparison with lattice data (section lv| and to study its behavior as the system approaches chiral restoration. 
In Appendix [A] we collect some of the SU(2) results, while in Appendix [b] we provide the detailed expressions for the 
thermal quark condensates to NNLO in ChPT. 

II. QCD CONDENSATES, SUSCEPTIBILITIES AND FOUR-QUARK CORRELATORS IN CHIRAL 

PERTURBATION THEORY 

A. General definitions 

Let us start from the QCD Euclidean Lagrangian including scalar sources: 

^QCD[q,q,s{x)] ^ q{i P - s{x))q-\ , (1) 

where the rest of the Lagrangian indicated by dots is irrelevant for our purposes, the (— , — , — , — ) metric is used and 
a sum over Nf flavor, Nc colors and Dirac indices is implicit. 
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The physical QCD Lagrangian corresponds to setting s{x) = Ai, the quark mass matrix. In the three flavor case 
A^=diag(mii,TOd,TOs), where rriu, md and correspond to the up, down and strange mass respectively. For simplicity, 
we will work in the isospin limit, so that to„ = — m and s{x) =diag(so(a;), So{x), Ss{x)). 

We will follow the external source method 16 to deal with the different two-quark (qq) and four-quark (qqqq) 
correlators of interest. Consider first the quark condensates at finite temperature: 



{qq)T = {uu + dd + ss)t ^ ^ tt-tt / • • • 99 exp / d'^xCQCD[q,q, s{x), ■ ■ ■], (2) 

'^QCDiy^l J Je 

where .^qcd[A^] is the partition function. In the above equation J^d'^x — Jq dr j d^x is the Euclidean (imaginary- 
time t = —it) version of the Minkowski volume i J d'^x, the averaging is performed over the thermal ensemble (an 
asymmetric box with imaginary time extension of /3 = 1/T ^ L and V — L^). The dots indicate the dependence on 
the rest of the QCD Lagrangian fields, not relevant for our purposes. Similar equations hold for the light condensate 
((9'?)i)T = {uu + dd)T- Recall that the light sector {u,d) is the most relevant one concerning chiral symmetry 
restoration, mostly due to the heavier strange mass. 

We consider the effective low-temperature representation given by Chiral Perturbation Theory (ChPT) [I4l - [l6] of 
the QCD generating functional, built from chiral symmetry invariance as an expansion in external momenta and 
quark masses: 

Zqcd[s\ ^ Z,ff[s]= Jvq^-e^pJ^d^xC,ff[r,six)l 

^eff = -^2 + ^4 + ^6 • ■ • , (3) 

We thus have: 



QGD\ 



( ^ 



^eff[M] \Sso{x) Sss{x) 



Zeff[s] 



(4) 



s=M 



where Ceff is the most general one made out of pion, kaon and eta fields </)", that respects the QCD chiral symmetry 
breaking pattern. These particles are the QCD low-energy degrees of freedom since they are Nambu-Goldstone bosons 
(NGB) of the QCD spontaneous chiral symmetry breaking. The subscript in the effective Lagrangian indicates the 
order in the ChPT derivative and mass expansion £ = 0{p'^*') over a typical scale ^ 1 GeV, and (jf" denote the NGB 
fields. Since the u, d, s quark masses are small compared with A^^^, they are introduced as perturbations, giving rise to 
the TT, K and 77 masses, counted as 0{p^). At each order, C^ff is the sum of all terms compatible with the symmetries, 
multiplied by chiral parameters, which absorb loop divergences order by order, yielding finite results. ChPT is thus 
the quantum effective field theory of QCD, and it allows for a systematic and model independent analysis of low- 
energy mesonic processes. The NGB fields are usually collected in the SU{Nf) matrix U = exp[iXa(f>"' / F], where, in 
the Nf = 3 case, Aq are the Gell-Mann matrices. The Lagrangian £2 is the non-linear sigma model: 



Co 



Tr[d^U''d''U + xiU + U'')] 



(5) 



with X = '2'Bqs{x), while F is the pion decay constant in the chiral limit. When s{x) — M, the following lowest order 
SU(3) relations hold: (qq) = -3BoF^, = 2Bom, M^j^ = Ba{m + m,) and M^^ = |Bo(m + 2™^). The ChPT 
power counting can be formally traced in terms of the counting in l/F"^ and so we will do in the following. The 
Lagrangians £4 and £§ are given in jl6j and [35, respectively, and contain the so-called low energy constants (LEG), 
Li and Hi (the latter are contact terms without NGB fields) for £4 and Ci for £g. 

The quark condensates (light and strange) can also be defined in terms of the free energy density z, which in the 
thermodynamic limit is given by: 



z = — lim log Z , 



(6) 



so that: 



m)i 



dm' 



(SS T 



dz 

diris 



{m)T = {{qq) 



IT 



iSS T- 



(7) 



4 



We turn now to the susceptibilities and their relation to four-quark correlators, which will play an important role 
in this work. Susceptibilities are defined as variations of the condensates with respect to the quark masses and 
are directly related to the thermal averages of four-quark operators measuring condensate fluctuations. Thus, the 
euclidean light scalar (or chiral) susceptibility is given by: 



(8) 



which relates the light susceptibility (the most relevant one regarding chiral restoration) with the four-quark correlator 
of the light combination {qq)i ^ uu + dd: 



Z[M] Ssoix) (5so(0) 



Z\s] 



(9) 



where T denotes euclidean time ordering, and so on for the strange, light-strange and full SU (3) susceptibilities: 



XAT) = -g^(^"«)T = = l^d^x [{T{ss)ix)iss)mT - {ss^t], (10) 

XUT) = -j^{{qq)i)T^-^{ss)T = -j^^^ [ d^X mqqUx)issm)T - {ss)t], (H) 

oiris am amoins 

X(T) - xiiT) + 2xisiT) + XsiT)= [ d^x [(r(gg)(a:)(gg)(0))T - (gg)?^] . (12) 

Note that, since the low-T representation of the free energy density is finite and independent of the low-energy 
renormalization scale /i [TT so are the quark condensates and susceptibilities, which can be expressed as mass deriva- 
tives of z. However, that is not the case of the four-quark condensates, which we define, for the different combinations 
of quark flavors, as: 



{{qq)a{qq)fi)T = lim(T(gg)a(x)(gg)^j(0))T, (13) 

where {qq)a stands for either {qq)i or (ss). Actually, at T = the low-energy representation of the four-quark 
condensates are divergent and scale-dependent, which is consistent with Renormalization Group (RG) analysis and 
holds also for other definitions of the four-quark condensates from the four-quark correlators, different from that in 



Eq. (13), as we showed in [5T] . 

Our previous discussion will allow us to relate, in the finite-temperature case, susceptibilities to four-quark corre- 
lators, which in the next section will be calculated within ChPT. Susceptibilities can be calculated either directly as 
mass derivatives of the two-quark condensates, or from the four-quark correlators, as in Eq. (|8|, always yielding a 
finite and scale-independent result. Even though four quark-correlators are not strictly needed for the calculation of 
susceptibilities, we will nevertheless calculate them, since they are needed to test the factorization hypothesis in the 
thermal case, because we aim to determine their validity as order parameters, and because they allow for a direct 
consistency check of our susceptibility calculation. 



B. ChPT thermal four-quark scalar correlators and condensates. 

Here we will calculate the relevant four-quark scalar correlators, by taking the corresponding functional derivatives, 
as in Eq. ([9]), from the effective Lagrangian in Eq. ([3]) at a given order in ChPT. We also calculate the two-quark 
condensate Eq. Q to the same order, to check the factorization hypothesis. We have: 

\ 5so{x) Sss{x) I J, 



(14) 
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(m)(x)(g-g)(0)), 



r 



+ r 



Sso{x) 



5ss{x) 



5Ceff[s] SCeff[s 



\ Ssq{x) 6ss{x) 



s=M 



s=M 



(15) 



All our results can be expressed in terms of the leading order (free) thermal meson propagators Gf{x), with i — tt, 
K, rj. Using standard finite-temperature methods, one can separate the T = part. The divergent contribution is 
contained in the x — and T — part. We follow the notation of fTT for thermal functions: 



Gf(0) = 
9i{M,T) - 



G.(0) 
1 

2^ 



.2 1 



dp 



Ep ef^i^p 



1' 



(16) 



with Ep — ^Jp' 



AP and G,(0) = Mq,^^^ [1 - D/2] /(47r)^/^ the T = divergent part in the dimensional regular- 
ization scheme, which we will use throughout this work. The renormalization of the LEG and the quark condensates 
up to NNLO in ChPT are the same as at T = and are discussed in detail in |31j . 

and 



Now, from Eqs. (1141) and (15), using the Lagrangians in 



{T{qq){x){qq){0))T, nlo = {m)T,NLO^ 
{T{qq){x){qq){0))T., nnlo = {qq)T,NNLO 



], we obtain the following results: 



24(12^6 + 2L8 + H2)S{T)d^^-^\x) + K'^{x) 



(17) 
(18) 



where the NLO and NNLO are 0{F'^) and 0{F'^) + 0{F°) in the ChPT counting respectively and K'^{x) is the 
connected part of the four-meson correlator at leading order and finite temperature: 

K^ix) = {Tr{x)M^)cl)\0)M0))T,LO ~ {rrmamT,LO = 2 {SGlixf + AGlixf + G^(x)2) . (19) 

We have expressed our results as a function of the square of the thermal quark condensate {qq)Ti whose explicit 
expressions are given in Appendix [b| Eqs. (B4)-(B7). As it happened in the T — case ^31] j up to NLO the four- 
quark correlator is just equal to the square of the quark condensate, but to NNLO the connected one-loop contribution 
breaks such equality. For a detailed diagrammatic description of the different terms contributing to the four-quark 
correlator, we refer to [31], since the diagrams are the same at T ^ 0. 

At this point, it is important to remark that we provide the full ChPT SU{3) result for the two-quark condensates 
to NNLO including all the relevant meson interactions. Previous works at finite temperature only included the kaon 
and eta (and other massive states such as nucleons) as free fields in the partition function [17j . which is reasonable 
for heavier particles when T « M}^,^, since the heavy states are Boltzmann suppressed. Therefore, apart for the 
study of factorization, our complete NNLO calculation of the condensate will allow to test the effect of including in 
the diagrams contact interactions containing strange particles. 

Similarly to the previous analysis, we calculate separately the light, strange and mixed four-quark correlators, which 
also factorize up to NLO in the product of the two-quark condensates, whereas up to NNLO we get: 



{T{qq)i{x){qq)imT ^ {{qq) 



lIT 



Bi 



16(8^6 + 2L8 + H2)5{t)5^°-^\x) 
+6G^(x)2 + 2G^(x)2-f ^G^(x) 



O 



F2 



{T{-ss){x){s.sm)T 
(r(gg)z(x)(ss)(0))T = {{qq)M^s)T + B^, 



8(4L6 + 2Ls + H2)5{t)5^''-^\x) + 2Gl{xf + -Gj^{xf 



O 



F2 



6ALeSiT)S^''-'\x) + 2G\{xf + -G^(a;)= 



F2 



(20) 

(21) 
(22) 



C. Factorization breaking at finite temperature 



As discussed in the introduction, scalar condensates play a relevant role in QCD, since they are directly related 
to vacuum properties. Attending only to their symmetry transformation properties, quark condensates of arbitrary 



6 



order ((gg)") should behave similarly as the two-quark condensate under chiral restoration, since they transform as 
isoscalars and are built out of chiral non-invariant operators with the vacuum quantum numbers. Actually, it is 
well known that such quark condensates appear directly in QCD sum rules, through the Operator Product Expansion 
(OPE) approach |36| . In that framework, the following hypothesis of factorization or vacuum saturation is customarily 
made: 

m')-{i-Jj^)m\ (23) 

and similarly for other condensates, which in the large-A^c limit simply reduces to {{qq)^) = (qq)"^. The second term 



between brackets in Eq. ( 23 1 comes from the exchange of indices (including color) between the first and second qq 
operators. 

The use of the factorization hypothesis has been a much debated tool in order to estimate the size of higher order 
condensates in the OPE. Let us remark that in |31j we have shown that factorization of the four-quark condensate 
does not hold within the model-independent QCD low-energy regime provided by ChPT, at T = 0. Actually, once the 
two-quark condensate is renormalized at a given order, factorization breaking terms are divergent and dependent on 
the low-energy scale. Therefore, the scalar four-quark condensate is not even a low-energy observable at T = 0. This 
factorization breaking is consistent with previous RG analysis [37] . Nevertheless, in [31) we found that factorization 
holds formally if the large Nc limit is taken before renormalization, since factorization breaking terms are 0{N^'^) 
suppressed. 

In this paper we extend the analysis performed in ^31, to the T ^ case, since, among other reasons, this can shed 
light on the use of the four-quark condensate ((g(?)^)r as an order parameter of the chiral transition. In particular, 
in case factorization holds, the four-quark condensate should behave as an order parameter, melting at the same 
critical temperature as the two-quark condensate. We obtain readily the factorization properties of scalar four-quark 



condensates by setting x — > in our results in section II B| It is clear then that, at finite temperature, factorization 



does not hold either. In fact, since S^^ ^\(^) vanishes identically in dimensional regularization [38], we get from 
Eqs.([l3]) and 

^ 1 + {SGliOr + ^GliQf + G^(0)^) + 0{l/F% (24) 

{qq)T 

We can calculate again the light, strange and mixed four-quark cases separately in S'J7(3), namely: 

^^i^ = 1+2^ [iGliQf + GliOr + gG^(O)^) + 0{1/F% (25) 



pi 



GliOf + lG^iOr'^+Oil/F'), (26) 



{{11)1 iss))T _ , 1 f^Ttn\2 , 2 



{{qq)i)T{ss)T F^ 



(G^(0)'+gG^(0)2)+O(l/F6), (27) 



In view of these results, several remarks are in order: First, the factorization breaking terms at T ^ are divergent 
and independent of the LEG, as for T = [3T], which was expected since the finite temperature result for the four- 
quark condensates to this order is just obtained by replacing Gi{0) — ?> Gf{Q). Since the quark NNLO condensates are 
rendered finite by the renormalization of the 0{p^) and 0{p^) LEG [31j, this means that the four-condensate to this 
order is divergent, also at T 7^ 0. Even subtracting the T = {) factorization breaking term, the result is still divergent. 



since, from Eq. (16), there are terms in Eq. (24) proportional to Gi{Qi)gi{Mi^T). Recall that, as mentioned at the 
beginning of this section, the four-quark condensates cannot be expressed only in terms of mass derivatives of the free 
energy density, unlike the quark condensates. 

Second, in the chiral limit, and contrary to the zero temperature case, for T 7^ the factorization breaking terms 
GJ (0) in Eq. ( 24 ) are finite and do not vanish. Let us recall that in dimensional regularization the T = propagators 



Gi(0) = when Mi — 0, since they are proportional to Mf . However, the thermal part gi{Mi,T) is finite and 
non vanishing in the chiral limit. Thus, the chiral limit is the only case for which the four-quark condensate can be 
considered an order parameter. Actually, we have checked that {{qq)'^)T/ {{<i<l)f)o follows the same temperature melting 
behavior as {{qq)i)T / {{q<i)i)o for the temperatures where ChPT is reliable. Note that chiral symmetry restoration 
takes place formally only in the chiral limit. 



7 



Third, for T ^ factorization holds formally in the iVc — oo limit, as it happened for the T = case f3T], since 
= 0{Nc) and therefore factorization breaking terms are once again 0{N~'^) suppressed. 

Partial results on the non-factorization of the thermal four quark condensates in certain approximations also exist 
in the literature. For instance, using the soft-pion and chiral limits it has been shown in |39| that if one assumes 
factorization at zero temperature, it will be spoiled by the lowest order thermal corrections. Our results support those 
in |39j since, on the one hand, we have proved factorization for T = in the chiral limit, and, on the other hand, we 
have also found that factorization breaks due to T 7^ contributions. 

In addition, an analysis of medium effects indicates that factorization is expected to be broken by suppressed 
contributions due to particles with the same quantum numbers as the operator under consideration j^Hj- Hence, we 
formally agree with [40j , since the leading term in 1 /Nc for the scalar four-quark condensate in the pure pion gas should 
factorize due to the absence of crossed terms {TT\qq\0). Nevertheless, by dropping these subleading breaking 

terms, some combinations of four-quark condensates, which include that in Eq.(25l above, have been proposed in |41| 
as order parameters. However, as we have just commented, even if these divergences are formally l/Nc suppressed, 
the four-quark condensates in Eq.(25 1 are actually divergent. Therefore, for finite Nc, the use of the order parameters 
proposed in '41] would require checking the cancellation of the divergences in the combinations suggested in 41]. 

Thus, one of the main conclusions of the present work is that the four-quark condensates in Eqs.(24) to (27), cannot 
be used as order parameters of the chiral phase transition, because they are divergent. Their use in previous works 
is valid only formally for Nc — > 00 or the chiral limit, but not in the physical case. We point out that a breaking of 
four-quark condensate factorization at finite temperature has been reported also in '42' within the framework of QCD 
sum rules. 



D. Scalar susceptibilities within ChPT 



We now turn to the ChPT evaluation of the different susceptibilities defined in section II A As detailed in Eq. i8 



and Eqs. ( 10 ) to ( 12 ), the scalar susceptibilities can be calculated either as a mass derivative of a quark condensate or 



from the corresponding four-quark correlators. Since we have available both the condensates (displayed in Appendix 
Ib| and the four-quark correlators given in section II B we can obtain the ChPT susceptibilities in both ways, thus 



checking our results for the four-quark correlators. Recall that precisely the additional factorization-breaking like 
terms in the NNLO correlators in Eqs. (18) and (20)-(22| give nonzero susceptibilities of chiral order 0(1). That 
order correspond to the derivatives of the NLO condensates in Appendix [b] since the 0{F^) is mass independent. 
If we calculate the susceptibilities from the four-quark correlators, we check again that the 0{F'^) vanishes (NLO 
correlators) while from the NNLO correlators, taking into account that 



d^x [Gfix)Y 



dM? 



Gf(0), 



(28) 



we readily obtain the different ChPT scalar susceptibilities to 0(1) in the chiral power counting: 



X(T) = B^a[24{12Llip) + 2Ll{^i)+HU^^))-Uly^-^6lyK-^^^r, 
652 (M., r) + 832 (Mk, T) + 2.g2(M^, T)] +o' ^ 



XsiT) 
Xis{T) 



XiiT) = Bl 



F2 



16 (SLlifx) + 2Ll{p) + i/J(/i)) - 12^^. - Ai^K -^i^n 



6g2{AU,T) + 2g2{MK,T) + -.g2(M,„ T) 



= B' 



O 



16 



8 {ALli^^) + 2Ll{^i) + - Avk - -v,, + 2g2{MK,T) + -g2(M„r) 



O 



p2 



MLliti) - AvK - -vr, + 2g2{MK,T) + -g2{M^,T) 



O 



with: 



1 



1 



92iM,T) 



327r2 
dgi{M,T) 
dM^ 



47r2 



dp 



Ep ef^^p - 1 ■ 



(29) 

(30) 
(31) 
(32) 

(33) 
(34) 
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A further check is that the results for the different susceptibihties are finite and independent of the low-energy scale 
/i, unlike the four-quark condensates, with the same renormalization of the low-energy constants, provided in [31j . 
which ensures that quark condensates are also finite and scale independent to that order. 

Note that, in the T ^ Mk regime, the thermal pion loops, i.e. 32 (-^-n- , 2^) , dominate over other particle thermal 
loop contributions, which are Boltzmann suppressed. An even more interesting regime is AIk ^ T 3> Af^, because 
it is related to the critical behavior [43 . Within our approach we cannot reach high temperatures but we can study 
the near chiral limit, where these functions behave as g2{M^^T) ~ T / {SttM^). We will remark this expected linear 
behavior when we plot our results below. 

Note that pion terms show up in the light susceptibility xi but not in the strange and mixed ones Xs , Xis which are 
therefore subdominant compared to the light one at temperatures below the transition. We have explicitly checked 
that in the limit — > 00 (A/^.r; ^ 00), the SU{3) light susceptibility reduces to the SU{2) result in Eq. (A5| in 
Appendix [Aj once the identification between SU{2) l\ and L\ LEG is made, as given in [TO] . 

The thermal scaling of the light susceptibility near the critical point {mu,d, T) — >■ (0"*", T^) reveals important features 
about the nature of the phase transition and is the subject of detailed analysis in lattice simulations [7]. The leading 
pion mass dependence from Eq.(30l near the chiral limit is Xit=o/-^o ~ "sf^ 1os(-^^/a*^) a-ud {xi,T — Xlt=o)^ ^/ Bq ^ 
^^jj- for ^ T ^ MK,-q- This leading behavior was already found in [35]. In addition, our results for the light 



susceptibility are consistent with a recent and model-independent ChPT analysis [33] of the mass, temperature and 
flavor dependence of the light susceptibilities, where a separate analysis of the quark connected and disconnected 
contributions was provided. 

After providing the analytic expressions for the susceptibilities, and addressing some formal aspects, let us discuss 
their phenomenology in connection with that of the quark condensate. 



III. PHENOMENOLOGICAL RESULTS IN CHPT 



A. Higher order meson interactions and the quark condensate 

As we have stated in the previous section, we provide in this work for the first time the full NNLO ChPT SU(3) 
{qq)T results, which are model-independent and given in Appendix |b] (the SU(2) one is also given in that Appendix). 
This result includes all the meson-meson interactions for tt, K, rj up to that order. At NLO, only tadpole contributions 
proportional to gi{Mi,T) appear, which is equivalent to considering an ideal gas made of the mesonic components. 
Previous studies of the condensate [T7] actually considered the contribution of the heavier K, -q states as free, while 
keeping higher orders in pure pion interactions. Here we will extrapolate our results to the critical point and make an 
estimate of the effect of considering interactions involving those strange degrees of freedom, by comparing with the 
NLO and with the pure SU (2) result in which kaons and eta have been decoupled. 

In order to study the thermal effects it is customary to normalize the quark condensate to its T = value and so we 
will do in what follows. An additional advantage of this normalization is that the ratio is QCD Renormalization Group 
(RG) independent, because the Bq global factors cancel. Let us note that, to NLO no LEG appear in (??)t/(9'?)o- To 
NNLO, only the 0{p^) LEG do appear. For the GhPT plots that we will present next, we will use the SU{i) set of 
L[ LEG from |45. and their error bands. An important remark is that, by definition, the H2 constant cannot be fixed 
by experimental meson data. This constant appears at NNLO and we estimate it as H2 = 2Lg as suggested in 
using scalar resonance saturation. In addition, as it is customary, at each given order we use F = F.^{1^ 0{Mf / F'^))^ 
= Mf{l + 0{Mf / F"^)) [16 taking the physical values for F.^ and Mf and including the corrections in the next 
order. 

Now, to comment on the relative size of different effects under discussion, we are going to extrapolate our results 
up to temperatures beyond the strict applicability limit of GhPT, typically estimated around T — 150 MeV p^. The 
reason is that it is much easier to explain the changes due to different effects in the curves of the condensate ratios, 
by comparing the points where their corresponding extrapolation vanish. We will refer to this point as the "critical 
temperature", T^, in the clear understanding that this is just for the sake of comparison between curves, since the 
particular value of this "extrapolated" temperature is just a very crude extrapolation and only the low temperature 



part is reliable and model independent. We will proceed similarly in section IV about the virial expansion. 

In order to study the restoration of the spontaneously broken chiral symmetry it is customary to use the non- 
strange condensate. The reason not to include the thermal strange quark condensate is that it has a very slow 
decrease as T increases, because its evolution is dominated by the larger mass, which is an explicit and not a 
spontaneous breaking. Thus, in the upper panel of Fig. [T] we show the thermal dependence of the normalized non- 
strange condensate within NLO and NNLO SU{2) GhPT, which can be compared with the SU{3) version in the lower 
panel. The gray bands surrounding the NNLO calculation cover the uncertainties in the LEG. We see that the NNLO 
correction is relatively small compared to the NLO one below Tc- The approximate temperature at which the SU(2) 
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FIG. 1: Non-strange quark condensate in ChPT as a function of temperature in the different cases explained in the main text. 



condensate vanishes is consistent with |17j , when we compare to the same order of approximation, taking into account 
that we use a more recent set of LECs. The reduction in Tc from the SU{2) to the SU{3) case can be interpreted 
as a "paramagnetic" effect coming from the increase in entropy (and therefore disorder) due to the addition of the 
strange degrees of freedom to the system. To NLO, the reduction on Tc from SU{2) to SU{3) is of 30 MeV whereas 
at NNLO is of 28 ± 12 MeV. This is the behavior expected when adding more degrees of freedom as they become 
more relevant near the phase transition, and it actually provides a natural explanation to the smaller values of Tc 
obtained in approaches such as the HRG which take into account all the relevant (free) degrees of freedom. In the 
SU{2) analysis we also observe a change in the value of Tc due to considering interacting pions (i.e., from NLO to 
NNLO) of about AT^ ~ —10 MeV. However, the result is inconclusive for the SU{3) case due to the uncertainties 
from the LEG dependence at NNLO. 



B. ChPT susceptibility 

We can also analyze the thermal evolution towards chiral restoration by studying the chiral susceptibility thermal 
dependence. In Fig. [2] we plot the ChPT xi,T result obtained for the SU{2) and SU{3) cases using Eqs. (301 and (A5|, 



respectively, normalized to their T — values and for the same set of parameters than the quark condensate in Fig. [Ij 
Note that we plot the ratio of the thermal to T = susceptibilities, which, once again, is QCD RG independent due 
to the cancellation of the Bq factors. We do not expect our low-energy analysis to reproduce the dramatic growth just 
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FIG. 2: Non-strange scalar susceptibility as a function of temperature for the SU{2) and SUiS) cases in ChPT at NNLO. The 
uncertainty bands, due to the in the LEG, have similar size for cases cases. The curves for the virial approach give a crude 
estimate of several sources of systematic uncertainties discussed in the text. 

below the critical point for xi- As it also happens for the quark condensate, we expect only to reproduce reasonably 
the low T side of the critical curves. We will provide a quantitative comparison with lattice values in section fV] below. 
As we have discussed above, the temperature evolution is governed by the T-increasing functions 32, so that, as we 
see m Fig.|2] the ChPT growth is roughly linear for temperatures high enough, according to the behavior in the chiral 
limit. In fact, although we have taken into account the full contribution of i:,K,r] loops in the SU{3) massive case, 
the kaon and eta thermal contributions provide a weak dependence, so that the SU{2) and SU{3) curves are very 
similar to one another up to the critical point, unlike the "paramagnetic" shift for the case of the quark condensate. 
In fact, although in both cases the SU (2) and SU (3) differ by terms of order e"*^^ ''/-^, the thermal dependence with 
T = T/Mk.-q is much softer in the susceptibility, since it comes from g2{MK.ri,T) ^ ^fre^^l'^ , than in the condensate 
where gi{AlK.ri,T) ^ r^(l/-yr)e"^/'^ for r <C 1. Again, the strange quark susceptibility is not showed, Xs/xo being 
also increasing but remaining very close to unity for the range of temperatures showed in Fig. [2] 

IV. THE VIRIAL APPROACH 

In previous sections we have been able to estimate the effect of the uncertainties in the LEG. However, ChPT is 
a perturbative expansion that, when truncated, neglects higher order corrections and cannot reproduce resonances. 
This is a first motivation to use the virial expansion, because it allows for a simple implementation of unitarized ChPT, 
which includes the numerically relevant higher order effects. An additional motivation to use the virial approach is 
that the interaction part can be clearly identified and described realistically. Both features will allow us to establish 
a consistent comparison with our previous standard ChPT approach. 

The virial expansion is a simple and successful approach already applied to describe many thermodynamical prop- 
erties of dilute gases made of interacting pions [551 [Ml HH] and other hadrons [T71 [2Z1 IMl HSl EHl HZ] • For most thermal 
observables, it is enough to know the T = scattering phase shifts of the particles that compose the gas. In principle 
these phase shifts could be taken from experiment avoiding any model dependence, so that it would not be necessary 
to go through the technicalities of finite temperature field theory. However, if one is interested in chiral symmetry 
restoration, and hence in scalar susceptibilities and quark condensates, one needs a model-independent theoretical 
description of the phase shifts in order to obtain their quark mass dependence, which cannot be obtained directly 
from experiment. Using directly one-loop ChPT for scattering lengths provides a remarkable description of the low 
energy hadronic interactions and should be accurate enough at very low temperatures [171. However, for temperatures 
further away from the threshold region, a more precise description of the scattering is needed, which in particular 
takes into account the loss of unitarity, and the absence of resonances, in the pure ChPT expansion. For that purpose, 
we will make use of the so called unitarized ChPT, at the expense of loosing the systematic ordering of the effective 
approach. 

The thermodynamics of a system of hadrons is encoded in the free energy density z = eQ — P, where eg = z"^^^ 
and P is the pressure. In the present work, we are interested in a multi-component interacting relativistic gas made 
of pions, kaons and etas in thermal and chemical equilibrium, so the pressure only depends on temperature T. Thus, 
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the second order relativistic virial expansion of the pressure reads [T7| 1211 [55] '■ 

z?^ = E ( ^^^^^ + ^1''^' + E ^'f ^'^^ I ' (35) 

where i = tt,K, j], the fugacities = e^*^*^* , Mi is the mass of the i species and the Bi and Bij are the virial 
coefficients for the gas. Expanding up to the second order in means that we only consider binary interactions. The 
coefficients 

/"OO 

(«) ^ ^ / ^pp2p-„/3(V5W-M,)^ (36) 

where the degeneracy is gi = 3, 4, 1 for tt, K, t] respectively, correspond simply to the virial expansion for a free gas 

dpp' log [l - e-^(V^)] . (37) 

The above free result is nothing but the HRG approach mentioned in the Introduction, when considering only the 
lightest hadrons tt, K, rj. Thus, the virial expansion provides naturally the corrections to the HRG due to interactions, 
which appear through the S'-matrix. For the meson-meson interactions, relevant for this work, this can be recast in 
terms of the elastic scattering phase shifts [T71 UHl [2H] • In this way, we can write 



nint ^« ^3 



^ / dEE^K,{E/T)A,,iE) (38) 



IMi + Mj 

where Ki is the first modified Bessel function and 

A^AE)=Y.{2I + l){2J+l)S'lJ^,{E), (39) 



LJ.S 



where the j § are the ij — ij elastic scattering phase shifts (chosen so that 5 — Q &t threshold Eth — Mi + Mj) 
of a state ij with quantum number /, J, S (isospin, angular momentum and strangeness), that we will explain below. 
The virial expansion breaks down typically where the dilute gas expansion does, for T ^ 200 — 250 MeV [2S]. In that 
regime 3> ^ C») so that the density of higher mass states an their interactions are Boltzmann suppressed with 
respect to pious. Hence, for our purposes of estimating systematic uncertainties, the ij = ttK and ttt] states can be 
neglected against the tttt interactions and we can drop the S index. 
From Eq. ([t]) we can then express the quark condensate as: 

m.)T^^^ = iomuo)-£-, (40) 

with a — l,s, mi = m and {0\{qq)a\0) — {{<j<l)a)T=o — deo/dnia. We emphasize again that in order to calculate 



Eq. (40) we need the dependence of S{E) on the quark masses as well as the vacuum expectation value. For that 



information we turn to ChPT in order to translate Eq. (401 in terms of physical meson masses: 



((.-.).)T = (0|(.-.).|0)(^1 + E^^). (41) 

with: 

cf = -F'g^{0\{qqUO)-\ (42) 

for which we will take the one-loop ChPT expressions from [121 130] • Since they depend on the ChPT LEC, we take 
the same values as in previous sections, obtaining; 

4'? = 1.02, eg = 0.59, c9« = 0.52. (43) 
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These numerical values are almost identical to those obtained in [50] , 

The behavior of the light and strange condensates within the virial expansion has been analyzed in detail for the 
meson gas in [^Hl I^MSO] and also including baryon interactions in [37| . In this paper we are interested mostly in 
four-quark condensates and susceptibilities. The four-quark condensates, e.g Eq. (|9|, cannot be obtained directly 
from the pressure as mass derivatives and hence the virial approach cannot give further information about them. The 
susceptibilities, defined in section [ll A| can be calculated taking one more mass derivative. Hence, 



dniadmp 

and again translating it in terms of meson masses: 



dmadmp ' 



1 



(0|(gg)„|0)(0|(g-g)^|0) 



q^dP_ 



(44) 



(45) 



By considering the ratios {{qq) a) t / {{<iq) a) o and Xa0{T)/xap{O), we cancel the overall Bq factors in the ChPT 
expressions for susceptibilities and condensates. These ratios are the quantities we will show in our plots. Note that, 
as it happened in the pure ChPT case, they are once again independent from the QCD renormalization group scale 
and can be expressed only in terms of meson parameters such as LEG, meson masses and decay constants. For these 
ratios we still need the values of {0\{qq)a\0) / Bq, for which we will take the one-loop ChPT expressions from [T^], and 



Xa/3{Q)/ Bq, which can be easily obtained from Eqs. (29 1 to (32) taking (?i = (?2 = 



Recall also that the contribution to the susceptibility of theiree part of the pressure Eq. ( 37 ) is precisely the same 



as the leading order ChPT results given in section II D so that the size of the interaction contribution from Eqs. (38) 



and ( 39 ) is a measure not only of the convergence of the virial series but also of the robustness of the pure ChPT 



contribution. 



Finally, in order to evaluate the interaction part of the virial coefhcients, Eqs. (38) and (39), we need the theoretical 



description of the meson-meson elastic scattering phase shifts, which are nothing but the complex phase of each 
scattering partial wave tjjs- These partial waves are obtained as the projection of the scattering amplitude in states 
of definite isospin /, angular momentum J and strangeness S. Let us remark that the unitarity of the S matrix 
implies that, for physical values of CM energy squared s, partial waves tjj for elastic meson- meson scattering should 
satisfy: 



Im tijs = cr\tijs\^ 



Im- 



1 



tiJS — 



1 



(46) 



where a = 2p/yG, and p is the CM momenta of the two mesons. Note that unitarity implies 



iz.5 = ^^e'^-, (47) 
a 

where Sjjs are the phase shifts needed in Eq. (39). The above equation leads to \tjj\ < 1/cr, and a strong interaction 
is characterized precisely by the saturation of this unitarity bound. 

In what follows we will first explain why the results for the susceptibility obtained combining standard ChPT with 
the virial expansion are very uncertain at very low temperatures due to a huge cancellation between the interactions 
in the scalar channels, and not reliable at moderate temperatures, since they do not describe data above typically 
500 MeV. Later on we will explain how this can be solved using unitarized elastic ChPT, which provides a realistic 
description of data up to roughly 1 GeV. 



A. ChPT Phase Shifts 



Let us then start discussing the partial waves tjjg, which arc obtained within standard ChPT as an expansion in 
even powers of momenta and meson masses. Dropping for simplicity the US indices, we find: t{s) — t2{s)+t4{s) + - ■ ■ 
where tn{s) = 0{p'^). Let us nevertheless recall that the ChPT series in only valid at low energies compared with 
47rF ~ 1.2 GeV, and in practice it is limited to scattering momenta of the order of 200-300 MeV above threshold, or 
about 400-500 MeV in energy. The reason is that, experimentally, for larger momenta several partial waves become 
resonant, a behavior that cannot be reproduced with a power expansion in energy. 
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The energy integrals that define the virial coefiicients, Eq. (38), extend to infinity, where the low energy ChPT 
expansion is no longer valid. In principle, one may think that this is not a severe problem because, for the temperatures 
we are interested in, the very high energy region should be suppressed by the thermal Bessel functions in the virial 



coefficients. However, as already explained, in the interaction part of the virial coefficients, Eq.(38), the tttt scattering 
phase shifts appear through the combination 

A^n{E) = 6ooiE) + 5S2oiE) + 9Sii{E) + ... (48) 

where we have omitted waves with J > 2 because they are significantly smaller than those with J < 1 below 1 GeV, 
which is the region of interest for our calculations. It is now important to remark that, as we show in the top left panel 
of Fig. [3] there is huge cancellation between 6oo and 5620- This cancellation was already observed for the scattering 
lengths in [T^. Using more recent determinations of the scattering lengths, we find aoo + 5a2o = 0.002 ± 0.009 versus 
aoo = 0.220±0.005 from gS], or 000 + 6020 = O.OlOiO.OlS versus ooo = 0.220±0.008 from 49J, namely, a cancellation 
of more than one, possibly two, orders of magnitude. Note that the latter pair of numbers are just a data analysis 
that does not use ChPT. 

As a consequence of that cancellation, the resulting A^^r is basically given by which at very low energies is 
extremely small. A similar cancellation also occurs both for the first and second mass derivatives of the phase shifts. 
Thus OAtttt /dAI-^ and d'^Aj^T^/dM^ are strongly dominated by the contribution from the vector channel, as we show 
in the left middle and left bottom panels of Figj3] Let us nevertheless remark that the J = 1 wave is suppressed at 
low energies by a q^'^ factor, and by itself yields a very small contribution to Atttt or its two first derivatives near 
threshold. Therefore, the model independent prediction of ChPT is that, at very low energies, the interaction part 
of the virial coefficient will be much smaller than naively expected from the size of each individual J = wave. In 
other words, the results of the virial approach with standard ChPT interactions at very low temperature follow the 
free gas approximation much closer than naively expected. 

Unfortunately, within ChPT we can state little more than the smallness of the interaction part of the virial coefficient 
at very low temperatures, because its value cannot be pinned down with numerical precision, since at very low energies 
what is left of A^^^ after the cancellation is even smaller than the size of higher order corrections. Still, in the next 
section, we will try to estimate those higher order contributions by means of unitarized ChPT. 

Nevertheless, one could think about providing ChPT results for moderate temperatures, say 100 or 150 MeV, where 
the bulk of the contribution to the integral extends beyond the threshold region and reaches, say, 500 MeV. However, 
in view of the left column in Figj3j it is clear that this can be done for the pressure and its first derivative, i.e. for 
the condensate, because A^^r and 9A^t/9M^ are growing functions of the energy, and soon enough the Su becomes 
dominant and the large systematic uncertainties near threshold become less relevant. But this is not the case for 
9^A^^/9M^, which is needed for the susceptibility calculation. Actually, as seen in the left bottom panel of Figjsj 
this second mass derivative is a decreasing function of the energy, so that the region nearby threshold, with its huge 
systematic uncertainties, continues to be dominant even for moderate temperatures, within ChPT calculations. 

Of course, if we consider even higher temperatures we find the usual caveats for a standard ChPT calculation, but 
this time even more severe. First of all because due to the cancellation commented above, the dominant contribution 
comes from the channel, which is also suppressed near threshold, and we can only find a sizable contribution when 
this 11 channel becomes sufficiently large, but that only happens around the p(770) resonance region. Unfortunately, 
whereas the ChPT description of the (/, J) — (0, 0) channel is fairly good, at least qualitatively, up to energies as high 
as 700 MeV, the presence of the p(770) resonance in the (J, J) = (1, 1) channel is not reproduced even qualitatively. 



Furthermore, the virial interaction coefficients in Eq. (381 are derived for an exactly unitary S matrix |25 , 26 , and 



this high energy thermal suppression quickly takes place for physical amplitudes respecting unitarity. But at this 
point we recall that ChPT scattering partial waves only satisfy unitarity perturbatively, i.e: Imt2 = 0, Imt4 = crt|,... 



and the unitarity constraint in Eq.(46l is badly violated pretty soon. Moreover, this violation grows fast with 
increasing momenta or in the vicinity of resonances. Hence, the contribution of ChPT is unphysically large due 
to the unitarity violation of the perturbative expansion and the contributions from the energy region E > 1 GeV, 
where any extrapolation is meaningless, become sizable already at temperatures of 150 MeV, for the ChPT virial 
calculation. In other words, for small temperatures, say formally T <^ the relevant momenta in the virial integrals 
are p '-^ \/Mt^T and hence the amplitudes are probed in the range p <^ M^^ where ChPT can be trusted. However, 
for T>Mt^ momenta are of the order of 2M7r or greater and the extrapolation of standard ChPT amplitudes does not 
describe meson scattering data. 

Within one-loop ChPT, the correct low energy expansion of the non-unitarized phase shift is Spfu ~ 17(^2 +Ret4), and 
gives a reasonable description of experimental data for very low energies with our choice of LEC. Then, for E ^ M^, 
Snu{E) ~ E"^, d6Nu{E)/dAI.^ ^ E"^, which produces additional powers of T in the interaction part compared to the 
free contribution, giving a large but unphysical weight to the higher energy contributions to the virial coefficients. 



Note in turn that by this same argument second derivative terms are subdominant in Eq. ( 45 ) with respect to the 



first derivatives. A rough estimation of the asymptotic behavior with T, say formally T ^ Af^, can be obtained by 
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FIG. 3: From top to bottom, the tttt scattering phase shifts {21 + l)Sij and their first and second mass derivatives compared to 
the combination A-mr = Soo + 5^20 + 9(5ii and its mass derivatives. On the left column we plot the ChPT results and on the right 
one the unitarized ChPT calculations. Note the huge cancellation that occurs, in all graphs and irrespective of unitarization, 
between the (/, J) = (0,0) and (2,0) contributions. Hence, A^^ and its derivatives are dominated by the (1,1) contribution, 
with only a few exceptions around the threshold region where the (1, 1) channel suffers an additional suppression. 



looking at the E ^ M^r behavior of the integrand in Eq. (38). The thermal function x'^Ki{x) weights the region 
a; ~ 1 so that asymptotically we can just trade factors of E in the phase shifts by T, which also allows to compare 
the interaction part in Eq. (38) with the free contributions Eq. (36). 

For all of these reasons, next we will make use of the so-called unitarized ChPT, which, under some reasonable 
approximations, combines ChPT at low energies with dispersion relations, and provides a realistic description of the 
data, reproducing the resonances relevant for this work, without introducing any spurious parameter other than those 
of ChPT. Of course, the price to pay is the loss of the systematic ChPT approach. 



15 



B. Unitarized interactions 



In order to check the influence in the virial integrals of higher ChPT orders, the violation of unitarity and the 
lack of resonances, we will extend the ChPT amplitudes by means of unitarization, and in particular, we will use the 
elastic Inverse Amplitude Method (lAM) [151 HHj which provides a remarkably good description of the meson-meson 
scattering data up to roughly 1 GeV. 



Unitarization methods provide amplitudes up to higher energies by using the fact implicit in Eq. (46), that the 
imaginary part of the inverse amplitude is known exactly. Naively, we can then impose the ChPT constraints to the 
real part of Rei""'^ ~ ^2^^(^2 + Rei4 + ...) to find that 

t = I ~ (49) 

Ret-i -ia 1 - t4/t2 ' 

This is the one-channel lAM [TBl - BD] . Although the use of the ChPT series in this naive derivation is only valid 
at low energies, the lAM can be derived also from a subtracted dispersion relation for the inverse amplitude, which 
justifies its extension to higher energies and even the complex plane. The details of the dispersive derivation can be 
found in [191 120) , but for our purposes here it is important to remark that the elastic cut is calculated exactly thanks 
to unitarity and the subtraction constants are calculated with ChPT, which is well justified since they correspond 
to evaluating the amplitude at very low energies. The lAM equation is valid at any energy as long as the left cut 
integral is well approximated by its low energy ChPT expansion, which is justified due to the subtractions, and as 
long as the energy where the amplitude is evaluated is sufficiently far from the inelastic region. Other terms due to 
so-called Adler zeros have also been shown explicitly to be negligible [5^ in this region. 

In summary, following a more rigorous derivation than the naive one above, the very same Eq. ( [49^ is recovered 
not only at low energies, but for most of the elastic region. Remarkably, this simple equation is able to describe 
meson-meson scattering data enlarging considerably the energy applicability range [19], while still reproducing the 
ChPT series at low energies. In addition, the lAM generates [T^ the poles in the second Riemann sheet associated 
with the resonances, from first principles like unitarity, analyticity and the QCD chiral symmetry breaking, without 
introducing these resonances by hand or any spurious parameter beyond the LEC of ChPT. Thus, the low-lying meson 
resonances are well described with this method, in very good agreement with the existing data. All these features 
can be reproduced with values of the ChPT parameters (LEC) that are fairly compatible with the values obtained 
within standard ChPT, despite being obtained from a fit to a much larger energy region. Nevertheless is important 
to point out that due to the nature of the lAM approach, the LEC needed to fit data and resonance poles with these 
unitarized amplitudes are only approximately those of ChPT. For this reason, we will use for the lAM phase shifts 
in the virial expansion the set of LEC obtained by fitting both scattering data and lattice results on meson masses, 
decay constants and scattering lengths. For the SU(3) case we take the values in [33] and for the SU(2) case from [51] . 
Note that we will only unitarize the pion-pion scattering amplitude since, as explained in previous sections and as it 
can be seen in our figures, this is the dominant contribution from meson interactions and, as we will see, correcting 
it with the lAM gives a considerably larger effect than including or not the kaon and eta interactions. Actually, and 
for the sake of simplicity, we will only consider two unitarized situations: either SU{2) or SU{S) but considering only 
free kaons and etas. With a lesser degree of rigor, the lAM can even be extended to the inelastic region ^50] above 1 
GeV, although that regime is not relevant for this work and the elastic formalism is enough for our purposes. Thus 
we prefer to rely on the most rigorous elastic formalism obtained from dispersion theory. 

Thus, in the right column of Fig. [3] we show the tttt phase shifts and their derivatives which are obtained from our 
unitarized ChPT (UChPT) calculations. Once again we find a huge cancellation between the (/, J) = (0, 0) and the 
(2, 0) contributions. In the case of the phase and its first derivative, which have less uncertainty in the cancellation 
and also become very small at threshold, the effect of this cancellation uncertainty is very small in the virial integrals. 
But this is not the case for d'^AT^^/dM^, as can be noticed when comparing the lower left and right panels, where 
we see that there are large uncertainties near threshold due to the higher order effects. Note that the change on each 
individual wave due to unitarization is rather small, but a mere 10% change in the (0, 0) channel produces a change 
of sign in the unitarized 9^A^^/9Af^. Hence, the interaction contribution to the susceptibilities is rather uncertain, 
but the overall uncertainty at very low energies is still small since the free contribution dominates by large. 

Note also that, for the unitarized case, we now draw the phases up to 1 GeV in order to show the almost complete 
dominance of the p(770) resonance contribution to A^^r and its derivatives above E = 500 MeV. The previously 
commented cancellation between the (0, 0) and (2, 0) contributions still persists at low energies, but deteriorates 
slightly above 500 MeV, where the p(770) contributions simply dominates because it is much larger than the others. 
This p{770) dominance is very relevant to asses the reliability of the UChPT results, since it has been recently shown 
that the p{770) mass dependence obtained with the one-loop lAM is in fairly good agreement with the most recent 
lattice calculations 15^. 
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FIG. 4: Non-strange quark condensate as a function of temperature in the virial approach, using Unitarized interactions in 
SU(2) (dotted Une)and in SU(3), only with pion-pion interaction and free kaons and etas (solid line) 



If now recall the result that in the narrow width approximation a resonance exchange contributes to the partition 
function as the free resonance state would do |25j, we conclude that the usual HRG with a free p{770), is, according to 
our results, a fairly consistent approach to include the tttt interactions. In contrast, one might naively expect that the 
a resonance, which is the nearest one to threshold and also has the quantum numbers of the vacuum, should provide 
the largest contribution to the susceptibility. However, and this is one of the remarkable results of this work, we have 
shown that in the threshold region it suffers a dramatic cancellation with the (/, J) = (2, 0) interaction. Therefore 
including just the cr as a free state in a HRG without the (2, 0) interaction, apart from ignoring the fact that the cr is 
by no means a narrow state, also neglects this very important cancellation. 

Moreover, the unitarized partial waves have a much softer behavior for large energies, namely t{E) behaves as a 
constant, giving rise to the asymptotic behavior 5u{E) ^ constant, d5u{E) / dM.^ ~ 1/E'^. Thus, for the susceptibility 
in Eq. ( 45 ) , the interaction part is suppressed with respect to the free one by inverse powers of T and the result is 
driven by the ChPT one. We have explicitly checked that the contributions to the integrals from energies higher than 
1 GeV are very suppressed now, and barely affect our results, contrary to the non-unitarized case. 

Thus, in Fig. [4] we plot the non-strange quark condensate using the virial approach with unitarized tttt interactions 
both within the SU(2) and SU(3) formalisms. The extrapolated melting temperatures are somewhat lower than those 
coming from standard NNLO ChPT calculations, already given in Fig.[l] particularly for the SU(2) case, this is partly 
explained, since as we have just discussed we are adding, in practice, the p(770) as an additional degree of freedom. 
Note also that the paramagnetic decrease between the SU(2) and SU(3) cases is just of the order of 6 MeV, which 
is smaller than the one estimated with NNLO ChPT. Nevertheless this smaller difference is less reliable since it is 
not calculated with SU(2) LEG obtained from those of SU(3), since the unitarized phases are obtained by fitting to 
different sets of data in both cases. 

Taking into account all the above considerations we have also plotted in Fig. [2] the chiral susceptibility with 
unitarized pion interactions. The results are similar to NNLO ChPT, and the difference provides the crude estimate 
of systematic uncertainties, giving rise to a quite consistent picture between ChPT and the virial approach. 



V. COMPARISON WITH LATTICE 



As explained in the introduction, the light scalar susceptibility is one of the main parameters analyzed in lattice 
simulations in order to identify its peak position as the transition point. A suitable quantity we can compare with 
is A(r) = m2 [xi{T) - xi{0)] /M^ = [xi{T) - x/(0)] /{iB^), which is given for instance in |1| for 2+1 SU{3) flavor 
simulations with almost physical quark masses. By subtracting the T = value, the lattice analysis of this quantity 
is free of ultraviolet divergences. Besides, A(T) obtained from the perturbative ChPT result in Eq. (301 is not only 
independent of Bq but also of the LEG (in particular of H2, which is subject to more uncertainty, as explained above). 
Thus, the ChPT result for A(T) depends only on meson masses and temperature. In the virial case, Eq. (45 1, there is 
no Bo dependence but the result still depends on the LEG, through the T = condensates, masses and phase shifts. 
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FIG. 5: Normalized relative non-strange scalar susceptibility in S(7(3) in terms of the relative temperature, in perturbative 
ChPT and in the unitarized virial approach. The points are lattice data taken from [6], where Tc ~ 155 MeV. 

As discussed previously, the lattice results predict a critical temperature considerably smaller than ChPT or virial 
extrapolations. This corresponds in part to the relevance of degrees of freedom of higher masses near T^. Qualitatively, 
one expects that those degrees of freedom produce a "paramagnetic" reduction of Tc due to the increase of entropy. 
Thus, in order to establish a more appropriate comparison with lattice results, we will represent the results in terms 
of the reduced temperature T /T^ for each approach (ChPT, virial and lattice) which is a way to compensate for the 
number of degrees of freedom involved. 

The results for A(T) are plotted in Fig. [5] At very low T, the curves remain close to one another. At high 
temperatures they grow roughly linear in T, up to the critical region. As stated before, this is a check of the 
robustness of both approaches, since the unitarized virial result is not only obtained within a different framework but 
it includes dependencies on the unitarization method and the LEC, thus giving an estimate of systematic uncertainties. 

Compared to the lattice data, we see that the two lowest points available are rather well described with the ChPT or 
unitarized virial curves. This is reassuring, since ChPT is meant to provide the low T model-independent tail. Once 
the temperature is re-scaled to Tc we see that the agreement is rather good even up to O.OTc, which is remarkable, given 
that ChPT captures only the qualitative features of the evolution towards chiral restoration but does not develop for 
instance a maximum, not even a sudden increase of A(r) near Tc. Actually, the lattice points reflect a clear departure 
from the ChPT prediction as they approach the critical point. 



VI. CONCLUSIONS 



We have analyzed several properties of the meson gas at low temperatures, regarding four-quark condensates 
and susceptibilities, within the ChPT and virial expansion approaches. Our analysis provides helpful results for 
understanding the behavior of the hadron gas formed after a relativistic heavy ion collision, below the chiral phase 
transition. 

The factorization hypothesis for four-quark correlators does not hold in ChPT at finite temperature to NNLO in the 
chiral expansion. This is an extension of a previous T — analysis and is a model independent result. In particular, 
it means that, in the physical case, the four-quark scalar condensate cannot be used as an order parameter for chiral 
restoration, since it contains divergent factorization breaking terms that cannot be renormalized. Nevertheless, there 
are two particular limits in which the four-quark condensate could be considered an order parameter, namely, the 
the large Nc limit, where factorization holds formally, and the chiral limit, where factorization is still broken, but, 
contrary to the T = case, just by a finite contribution. Let us remark that the factorization breaking terms are 
precisely those needed to provide finite and scale-independent scalar susceptibilities for the light and strange sector, 
including the mixed one. We have provided explicit expressions for those susceptibilities in ChPT to leading order in 
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their chiral expansion. The most important one, regarding chiral restoration, is the hght scalar susceptibihty, which 
grows hnearly in T at low temperatures. 

In order to establish properly the previous factorization results, we have calculated the two-quark condensate in 
5t/(3) up to NNLO in ChPT, including meson interactions with strange degrees of freedom {ttK, ttt], K-q and so 
on). In particular, this allowed us to discuss the effect of those interactions in the crude determination of from the 
extrapolated condensate. In addition, since the LEG of fourth order enter at that level, we have been able also to 
estimate the influence of the LEG in the GhPT determination of T^.. The net effect of strange interactions is about 
ATc = 28 ± 12 MeV where the error comes from the LEG uncertainty. 

An important part of our work has been devoted to the comparison of the GhPT approach with the virial or density 
expansion, especially for the light scalar susceptibility. We have shown that one should consider unitarized interactions 
when probing the low and moderate temperatures of interest for this work. Moreover, the unitarized virial curves 
remain close to the standard GhPT calculation, which is a reflection of the robustness of both approaches, at least at 
low and moderate T, and shows that the pion interactions are suppressed in thermal observables. This suppression 
is even larger than expected due to a huge cancellation between the scalar channels with isospin and 2, which is 
already observed in the data below 1 GeV, but we have found also to occur for the first and second mass derivatives 
of the interactions. The existence of this cancellation implies that if the sigma particle, which dominates the scalar 
isoscalar channel near threshold is included alone as a free state, as in the usual hadron resonance gas approach, it can 
give large deviations from the model independent GhPT approach even at very low temperatures. As a consequence 
of this cancellation, the p(770) dominates the interaction contribution starting at moderate temperatures of the order 
of 100 to 150 MeV, accelerating the melting of the condensate and the susceptibility growth. 

The comparison with lattice data shows a remarkable agreement at those temperatures, when they are re-scaled with 
respect to the corresponding critical temperatures, i.e., compensating by the number of degrees of freedom involved in 
the calculation. When we compare between the SU{2) and SU{3) cases, we do not see a significant difference between 
considering free or interacting kaons and eta compared to the effect of unitarizing the interactions. In the standard 
virial treatment, where the phase shifts are considered at T = and the thermal correction comes from the weight 
density functions, we have showed that the unitarized interactions decrease with inverse powers of T with respect to 
the free contribution. 



Appendix A: SU(2) results 



Here we collect, for completeness, the SU{2) results for the four-quark correlator, factorization and scalar suscep- 
tibility. 

For the four-quark correlator we get: 



{Tm{x){qq)mT , NLO ~ 



T, NNLO 



T, NNLO 



(Al) 
(A2) 



where we have defined Kgjj2{x) as the connected part of the four-pion correlator at leading order and finite temper- 
ature: 

K^u2i^) = {rr{^)u^)c^''mbmT.Lo - {rrmami.Lo = ^Gli^f- (as) 

As it happened in 5f7(3), to NLO the four-quark correlator can be expressed again as the square of the quark 
condensate, but not to NNLO where Kgjj2 breaks such factorization. As we did in the main text, we are giving 
Eqs. (Al) and (A2| simplified in terms of the explicit expression for {qq)T. nnlOj which are given in Appendix pi 



As for factorization, we have from Eq. ( |A2 

{im?)T 



^(G.(0)+5i(M^,r)f , 



(A4) 



which is again divergent with the standard GhPT renormalization 
The SU (2) susceptibihty is given by: 



(A5) 
(A6) 
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An important comment is that the renormahzed SU{2) LEG can be written in terms of the SU (3) ones by performing 
formally an expansion for large strange quark mass in a given observable calculated in SU{5) and comparing with the 
corresponding SU{2) expression |T1], although the numerical difference between them is small. In the SU{2) case, the 
only two combinations of LEG appearing in the NNLO expression for {qq)T/ {QQ)a are Zg and l"^ + h^, which can be 
readily expressed in terms of the SU{3) LEG. Applying such LEG conversion in Eq. (A6| one gets directly the SU{3) 
expression in Eq. (301 if the g{Mk^n,T) are neglected. 



Appendix B: Finite temperature quark condensates to NNLO in ChPT 



In this appendix we will provide the NNLO results for the two-quark condensates at finite temperature. As explained 
in the main text, the corresponding four-quark condensates cannot be obtained just by squaring these results, but one 
also has to add the factorization breaking contributions described in Eqs. (A2) and (18). The renormalization needed 



to render the quark condensates finite and scale independent is the same as for T — 0. Therefore, for all the technical 
aspects concerning the conventions for the needed £4 and LEG and their renormalization, we refer to [31j . 
For convenience we define: 



M2 M2 

0^ log °* 



327r2F2 



2p2 ' 



MT) = F 



327r 



1 ^+log.^'^"A 3-^^'-^^ 



(Bl) 



where we denote, following the notation in |16j . = /ii(0) = Mi (51 — 0). 

The final expressions for the two-quark condensates, finite and scale-independent, up to NNLO, at T 7^ are given 
by: 



,SU{2) 
l/NLO 



(T) = -2BaF^ i 1 



F2 



{hl + ll)-'ii,^{T)^, 
3 3 Af 2 



m)i)f^o\T) 



, -2B,F'l^l + ^[{H^^+AL\ 



1 

F2 



-iM^,^l^{T)vAT)+^- 



,,~ALl)Ml^ + 2{L\ 



24 

^fiAT) [(3LI + 2L;~6LI 

^^i^,{T) [{-iL\-2Ll 
2AMI 



-2Ll)M^^] 



+ 6Lg - 48L; - 12L^) + 2 (15L^ + 4Lg - 30Lg 
'5 ~ 2Lg — 2Z,g) Mq^ 



- 24L;) Ml 



(B2) 



(B3) 



.1 + 2LI) Ml + SLlMlj,] - 3m.(T) - 2^iK{T) - \h{T)^ , 



(B4) 



^2{Ll 



-2Lg)MoV] 



s 



7,)Ml + {2Ll 



-2Lg)MoV] 



' - -2{L\ ' 



8 (3L^ - 



'5 - m 



'6 + ^' 



- &LI - 48L; - 18i^) M^ + 2 {iFl - ALl 

'7 + -^s) 



- 6Lg + 48L! 



8)))M0J 



'qk\ 



^8 

+ ^ [(3^1' - 2(^2 + C^s) Mt + 4 (Cl - Cl) M^Ml 
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(B5) 



20 



{ss)nlo{T) 
{ss)nnlo{T) 



+ 



-BoF' { 1 

{ss)nlo 
32 



BoF^ 



^ [- {H; ~ ALl + 2Li) Ml + 2 (i/J + 4Lg + 2L^) M^^] - 4A.;,(r) - ^/i,(r)| , (B6) 



^J.K{T)^ln{T) 
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-M^^lJi^{T)vK{T) + - (16M2^ 
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■ 48L; + 24L^) MlM^K 



(B7) 



where the Gell-Mann-Okubo relation "iM^ = ^Mqj^ — Ml for the SU (3) leading order masses has been used and the 

renormalized L[, Z[ and c[, CI constants depend on the scale /i as explained in [3T]. We recall that L4 and L5 appear 
because of the meson wave function and mass renormalization. In S'J7(3), the constant Ly stems from the eta mass 
renormalization. 
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